Fidelity freeze for a random matrix model with off— diagonal perturbation 
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The concept of fidelity has been introduced to characterize the stability of a quantum-mechanical 
system against perturbations. The fidelity amplitude is defined as the overlap integral of a wave 
packet with itself after the development forth and back under the influence of two slightly different 
Hamiltonians. It was shown by Prosen and Znidaric in the linear-response approximation that the 
decay of the fidelity is frozen if the Hamiltonian of the perturbation contains off-diagonal elements 
only. In the present work the results of Prosen and Znidaric are extended by a supersymmetry 
calculation to arbitrary strengths of the perturbation for the case of an unperturbed Hamiltonian 
taken from the Gaussian orthogonal ensemble and a purely unitary antisymmetric perturbation. It 
is found that for the exact calculation the freeze of fidelity is only slightly reduced as compared to the 
linear-response approximation. This may have important consequences for the design of quantum 
computers. 
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I. INTRODUCTION 

The concept of fidelity was originally introduced by 
Peres to characterize the quantum-mechanical stabil- 
ity of a system against perturbations. Recently it enjoys 
a renewed popularity because of its obvious relevance for 
quantum computing. In the present context the works 
which focus on random matrix aspects are of particular 
relevance. First the paper by Gorin et al. has to be 
mentioned where the Gaussian average of the decay of the 
fidelity amplitude was calculated in linear response ap- 
proximation. For small perturbations the authors found 
a predominantly Gaussian decay, with a cross-over to ex- 
ponential decay for strong perturbations, in accordance 
with literature 0, Q . The results of the paper could be 
experimentally verified in an ultrasound experiment .Sj 
and in a microwave billiard 0. Using supersymmetry 
techniques, the limitations of the linear response approx- 
imation could be overcome, yielding analytic expressions 
for the decay of the fidelity amplitude for the Gaussian 
orthogonal (GOE) and Gaussian unitary (GUE) ensem- 
ble. Quite surprisingly a recovery of the fidelity was 
found at the Heisenberg time "t*] which was interpreted 
as a spectral analogue of an Debye- Waller factor |8|| . Ref- 
erence is the basis for the present work. 

The Gaussian decay observed for small perturbation 
is caused by the diagonal part of the perturbation in the 
eigenbasis of the unperturbed Hamiltonian. This was the 
motivation for Prosen and Znidaric to look for pertur- 
bations with zero diagonal, first in classically integrable 
systems 0- Later on they extended there studies to 
classically chaotic systems [l^. In linear response ap- 
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proximation they found an plateau in the decay of the 
fidelity. Only after extraordinarily long times the decay 
started again, exponentially below, and Gaussian beyond 
the Heisenberg time. It remained, however, an open ques- 
tion whether this freeze of the fidelity is reality or whether 
it is just an artifact of the approximation. It was the mo- 
tivation for the present work to answer this question by 
extending the previous supersymmetry calculation to the 
freeze situation. It will be shown that the freeze is present 
also in the exact calculation. This may have important 
consequences for quantum computing. If one succeeds in 
imbedding the atoms representing the qbits into an en- 
vironment coupled only via an off-diagonal perturbation 
to the atoms, an enhancement of the system's stability 
by orders of magnitude is expected. 

The present results are not restricted to random ma- 
trices. In reference it is shown that, e.g., kicked 
tops with a corresponding dynamics follow exactly the 
random matrix predictions of the present paper. 



II. THE LINEAR RESPONSE 
APPROXIMATION 

The fidelity amplitude is defined as the overlap integral 
of an initial wave function |^/;) with itself after the time 
evolution due to two slightly different Hamiltonians Ho 

and Hx = Ha + \V, 

A(r) = (^|e2-^^-e-2-««-|V.) , (1) 

where the time r is given in units of the Heisenberg time. 
For chaotic systems fx{T) is independent on the initial 
condition, and we may replace expression by its av- 
erage over 1^), i. e., 
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where N is the rank of the Hamiltonians, and the brack- 
ets denote an ensemble average. Under the assumptions 
that (i) Hq is taken either from the Gaussian orthogonal 
(GOE) or the Gaussian unitary ensemble (GUE) with a 
mean level spacing of one in the band centre and (ii) that 
the variances of the matrix elements of V are given by 
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Gorin et al. obtained for the fidelity amplitude in the 
linear response approximation 



/,(T)^l-eC(T), 
where e = Att'^X^, and C(t) is given by 
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^2, /?(■'■) is the two-point form factor, and P is the uni- 
versality index, i.e. (3=1 for the GOE, and /? = 2 for 
the GUE. For the Gaussian ensembles b2^{T) is known, 
and C(r) can be explicitely calculated j3- The range 
of validity of the linear response approximation can be 
somewhat extended, by exponentiating Eq. 
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The authors argued that the errors of the approximation 
should be fairly small for A ~ 0.1 and negligible for A ~ 
0.01 (corresponding to e = 0.4 and 0.004, respectively), 
which was fully confirmed by the exact calculations Q- 
The Gaussian decay for small perturbations is caused 
by the diagonal part of the perturbation. This is imme- 
diately evident from Eq. 10): For small perturbations V 
can be truncated to Vdiag, its diagonal part in the basis 
of eigenfunctions of i^o- In this regime Eq. ^ reduces 
to 
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where the £"„ are the eigenenergies of Hq 4] . This sug- 
gests to consider perturbations with vanishing diagonal 
matrix elements in the eigenbasis of Hq 12] . In the linear 
response approximation one then obtains 
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where Cfroozo differs from the expression ^ derived pre- 
viously only by the fact that the term r^//? is missing on 
the right hand side of the equation The resulting 

decay of the fidelity amplitude is extremely slow. It will 
be discussed below and compared with the exact result 
as obtained from the supersymmetry calculation. 



III. THE PURELY IMAGINARY 
ANTISYMMETRIC PERTURBATION 

To apply the supersymmetric technique of Reference 
the Hamiltonian needs to be invariant under the action 
of the orthogonal/unitary group. This is not the case 
for a GOE perturbation with a deleted diagonal. How- 
ever a purely imaginary antisymmetric matrix meets with 
both requirements, zero diagonal elements and orthogo- 
nal symmetry. Therefore it is an ideal candidate. We 
consider the Hamiltonian 



Hx = Ho+ iXV , 
where Hq is taken from the GOE, i. e. 



(9) 
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(SikSji + SiiSjk) , (10) 



and V is real antisymmetric, i. e. 

{VijVki)y = {6ikSji - SilSjk) ■ (11) 

The variance of the matrix elements has been chosen to 
have a mean level spacing of one for Hq and of order 
1/VN for V. Note that in contrast to Reference 7 the 
mean density of states does not remain constant with in- 
creasing perturbation, but decreases with increasing A. 
In fact it is irrelevant in the present context, whether 
a defolding to a constant mean density of states is per- 
formed or not. Such a defolding would imply an addi- 
tional factor of l/-\/l -I- {ttXY /N on the right hand side 
of Equation Q , which in the final limit N ^ oo reduces 
to one. These definitions are consistent with the normal- 
ization used by Gorin et al. 0- Expressing ftir), see 
Eq. |(2Jl, in terms of its Fourier transform 

Mt) = J dE, di?2e2-(^i-^^)- iEi,E2) , (12) 

we have 

R, (Ei,E2) oc — / tr ( 

^ ' ^ n\ \Ei^ - Ho - iXV E2+ - Hq 

(13) 

with E± = E ± IT]. We rewrite Re {Ei,E2) using the 
formula 
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d d Aet{A+J)dei{B + K) 



AB ^ dJnrndKmndet(A- J) det(B - K) 

J=K=0 

(14) 

In our case A = E2+ — Hq is real symmetric and B = 
— Hq — iXV is Hermitean. According to the univer- 
sality classes of A and B we write the determinants as 
Gaussian integrals over real and complex wave functions, 
respectively. We obtain 



d det(A + J) 



^ dJnm det{A - J) 

n rn. ^ ' 
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where the commuting integration variables are real. We 
adopt the usual convention and use latin letters for com- 
muting, and greek ones for anticommuting variables, re- 
spectively. For B we obtain instead 
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d det{B + K) 



dKmn det{B - K) 



= i / d[a]d[b]d[Tj]d[T]* 



K=t) 



(flmOn + bmhn - VniVn - VrnVn) 



(16) 



where Zi is complex, and and bi are its real and imag- 
inary part, respectively. Collecting the results we arrive 
at 
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^(x X +v V ~£ r)(a a +b b - v* n -n /e'^E„.,„ 

/ ^ \'^rf^m ^ ynyrn Sn'^^^i S^S.m 7 V^^^i"^ ^ ^m^n 'Ini In 'ira'in) \ ^ 



V 

(17) 



Ha 



r 



The commuting integration variables are all real. Now 
the average is taken over real symmetric Hq using 
Eq. (0 



■)ho ^exp 
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Str(LZ)^ 
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where Z is a supermatrix given by Z = 'Yin ^"^n with 
= {xn,yn,£.n,^n,an,bn,rin,ril), and L = diag(l4,-l4) in 
the advanced-retarded block notation. Z is exactly the 
matrix given in Table 4.1 of Reference lisl denoted by 
VWZ in the following. The matrix B has an orthosym- 
plectic symmetry, i. e. in Boson-Fcrmion block nota- 
tion the Boson-Boson block is real symmetric and the 
Fermion-Fermion block is Hermitean selHual. For the V 
average we obtain with Eq. ((111) 



exp (-A2Str(Kr)^ 



(19) 



Here T = Y,n ^"^n "^'^^^ = ^n, Vn.Vn) super- 
matrix K is given by 
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where we used the Pauli matrices 



(20) 
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In K the off-diagonal nature of the perturbation is en- 
coded. 

The subsequent steps are the same as described in 
0, ■ After transforming equations p8|) and (|19|l by 
means of two Hubbard-Stratonovich transformations, the 
integrations over the the a, 6, x, y variables, and over the 



auxiliary variables of one Hubbard-Stratonovich trans- 
formation can be performed resulting in 
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Here P = diag(l, 1, —1, —1). The matrix a has the same 
orthosymplectic symmetry as B and reads in advanced- 
retarded block notation 



<^RA <7RR. 



(23) 



Introducing the notation i?i/2 = EzLE/2, and substitut- 
ing aAA and a^R by aAA + E/2 and aRR ~ E/2, respec- 
tively, we obtain 



R,{E^,E2) cc 



4iV3 



d[<j]Sti (Pal^PaRA 

+ E Str{aAA-<yRR}] 



e 4JV I 

g^[Str(K<TRH)2-_BStrK(Tfl„] 

I -N/2 



[Sdct(cr - E)] " 



(24) 



This expression can be evaluated in the limit N oo 
by a saddle point approximation. 



IV. SADDLE POINT APPROXIMATION 

The next steps are a direct repetition of the corre- 
sponding ones in Reference Q We shall adopt the nota- 
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tion from VWZ, which is the main source for the follow- 
ing calculations. First we diagonalize cr, 



where i? is block diagonal and 



To 



VI + ^12^21 «il2 



-1^21 Vl + ^21^12 



(25) 



(26) 



(see VWZ, Eqs. (5.28+29)). The integration of the di- 
agonal variables of a can be performed by means of the 
saddle point approximation, au at the saddle point reads 



OD = 



SA 
SB. 



(27) 



where the advanced and retarded saddle points are given 
by 



SA/R = -{E±iA) 



and p is the density of states. In the following we shall 
restrict ourselves to the band centre, = 0, where Eq. 

reduces to SAjR = ±nr/N. We then have for the 
matrix a at the saddle point 



^ iV / l{l + 2ti2t2l) 2ti2Vl+t2ltl2 \ .„ql 
TT ^ 2^21^1 + il2t21 z(-l-2t21<12) ) ' ^ ^ 

The matrix R (see Eq. (|25|l ') does not enter, since it 
commutes with cr^i at the saddle point. We obtain 
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R,{Ei,E2) cx ^ hFiti2)d[ti2]Str{PaARP<JRA) 

g-lW [EStiiaAA-<rBB)-^^Sti{KaAAf] j-gg-j 

where the integral is over the elements of the matrix 
ti2, parametrizing the saddle-point manifold. The func- 
tion J^(ti2) — Sdet~^^^(l -I- ti2t2i) is the Berezinian of 
the coordinate transformation Eq. (|26|l . Using Eq. I|29() 
the various terms entering Eq. may be written as 



4iV 

StT{crAA- ctrr) = I Str(ii2i2i) 



(31) 



Str (Kcr_R_R) 



TT 

AT \ 2 

Str[K(l + 2ii2i2i)]'t32) 

2N\' 



Sti{aARP<JRAP) - ( — ) Str (i2iVl + ii2t2iP 



X tuVT+h^P) . (33) 

We proceed further by diagonalizing the matrices ti2 and 
t2i ■ This is achieved by the radial decomposition 



ti2 = U^'^MUU2 , t2i = U2^U''MUi 



(34) 



with diagonal M = diag(/ii, /i2, i/i, «/^). The matrix 
U = I2®?7isa4x4 block diagonal matrix, with U 
G SU(2) (see VWZ, Eq. (1.18)). The Ui {i = 1,2) may 
be parameterised as Ui — ViOi, where the Oi — Oi ® I2 
are 4x4 block diagonal, with Oi S SO (2). The parame- 
terisation of the Vi in terms of anti-commuting variables 
is postponed to App. ^ If we moreover introduce 



X = A/^ = diag(a::, y, —z, —z) , 



(35) 



we can write Eqs. (jSU to using Eq. (jSUl and Eq. 



m 

I StrX 

TT 

2 



(36) 



— Str[Ki(l + 2X)]\(37) 

TT 



2N 



Str ( ^/X^/lTXPl 



Str {(TAA - ctrr.) 
Str (KaRR)'^ 

StT{aARPa-RAP) 
where 

Ki = C/iKC/f 1 , Pi = [/iPC/f 1 . P2 ^ UU2PU2^U-^ . 

(39) 

Under the transformations Eq. (|34|1 and Eq. H35|) the 
measure transforms as 



XVl + XP^] , (38) 



d[ti2\ - g{X)dn{Ui)d^jL{U2)diJL{U)d[X] 



(40) 



The function Q has been calculated in VWZ (Eq. K.17). 
The average in Eq. H3U|I is over the elements of the ma- 
trices X, U, C/i , C/2 • Only P2 depends on the matrix ele- 
ments of 1/2- It will be shown in App. ^that U2PU2^ 
averaged over the matrix elements of U2 is nothing but 
a multiple of the four-dimensional unit matrix. Thus the 
U dependence cancels. We are then left with an aver- 
age over X, y, z, and the matrix elements of Ui. Inserting 
these results into Eq. H3UI) . we get 

R,{EuE2) cx ^ j T{X)g{X)d[X]dn{Ui) (41) 

Str [X{X + l)Pi] e-"SStrX-T^Str[(l+2X)K,]= ^ 

where we employed the definition e = 47r^A^, see above. 
The Berezinians J^{X) and Q{X) can be comprised in one 
measure function /i(A") which was given in VWZ: 



m(x) 



k - y\ 



z{l-z) 



^xy{x + l){y + l) [z - xf{z - yf 



(42) 



Substituting expression H41|l for Re{Ei,E2) into Eq. 
(P|l . and introducing E = {Ex ~ E2)/2 and E = 
[El -\- Ei)/2 as new integration variables, the E inte- 
gration generates a delta function, whereas the E inte- 
gration corresponds to an energy average. The result is 
(see Reference for details) 
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/^(t) oc — / du / / dzii{u,v, z)dij{Ui) 

^ Jo Jo Jo 
(5 (r - M - z) StT{X{X + l)Pi) 

exp(-^Str[(l + 2X)Ki]') , (43) 

with u = {x + y)/2. In addition we shall use = {x — y)/2 
as another new variable, and replace z hy t — u every- 
where, which is admissible because of the presence of 
the delta function. The integration domains of the ra- 
dial variables u, v and z are dictated by the hyperbolic 
symmetry of the saddle point manifold, i. e. we have 
non-compact integration domains < x,y < oo for the 
bosonic coordinates x, y and a compact integration do- 
main < z < 1 for the fermionic coordinate z. For more 
detail on this point, see VWZ. We still have to integrate 
over the matrix elements of t/i. 

V. INTEGRATION OVER THE GRASSMANN 
VARIABLES 

We recall that Ui — ViOi. Since Oi commutes with 
P and with K, the Oi integration is trivial, and we are 
left with the integration over Vi . The parametrization of 
Vi in terms of Grassmannian variables and the calcula- 
tion of the traces in equation 1)431) is quite involved, and 
is postponed to the appendices. Here we note only the 
results: 

Str[X(X + l)Pi] = 4:v{2u+l)B 

+2[2u{u t{2u + 1 - r) -f w^] 

+A[T{2u+l-T) + v\A-2a), (44) 

and 

^Str[(l + 2X)Ki]' = 
o 

I 



t{2u + 1 - t) - -I- 2{A - D){t^ - t;2),(45) 

where 

A^aa*+I3l3*, B ^ aa* - ^P* , D = i{aP* - (3a*) . 

(46) 

and a — aa*(3l3* . a, a*, /3, (3* are anticommuting 
variables. It follows 

g-^Str[(l+2X)Ki]" ^ ^-i[T{2u+l-r)-v^+2(A+D)(T'^-v')\ 

= [l-e{A~D){r^~v^)] 

These results are inserted into equation H43() . The mea- 
sure is given by 

d^i{Ui) = 27rd/i(yi) cx dada*df3d(3* . (48) 
Therefore only the terms proportional to a survive the 
integration over the antisymmetric variables. We obtain 

/e(T) (X ^ J ^,{u,v,z)SiT~u~z)e-i(^^'-+'-^^-^") 
[1 + e{T^ - v'^)][t{2u + 1 - t) + v^]dudvdz , (49) 

which is almost our final result. 



VI. RESULT AND DISCUSSION 

The final result is obtained by an VWZ-like integral 
(see VWZ, Eq. (8.10)) and is given in the present case 

by 



r u 

2 [ d [ ^'^^ (r — m)(1— t + m) 

J J\u-^ - vmu + - v^] («2 - t2)2 

Max(0,r-1) ^ 

X [1 + e{T^ - V^)][t{2u + 1 - t) + ^;2]g-§[r(2u+l-r)-«^] _ (gg) 



The constant of proportionality was fixed by the con- 
dition /e(0) = 1 (see Reference 8). The only difference 
of Eq. H5()(l to Reference 8, where a GOE perturbation 
was considered, is the additional factor [1 + e(r2 — v'^)] in 
the integrand, and a minus sign with the term in the 
exponent, where in the GOE case there is a plus sign. 



Figure n shows the fidelity decay for different pertur- 
bations, as calculated from Eq. H5U|) . together with the 
result from the exponentiated linear response approxi- 
mation. For comparison the fidelity decay for the case of 
a GOE perturbation '^) is shown as well. We see that 
the linear response approximation is able to describe the 
fidelity decay for quite a long time very well. For still 
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larger times the linear response approximation underes- 
timates the decay, as compared to the exact result, but 
still the decay is by orders of magnitudes slower as for a 
GOE perturbation. 

Figure [3 shows the fidelity (r) for three fixed val- 
ues 0.5, 1.0, 1.5 of T as a function of the perturbation 
e. The figure demonstrates that the freezing effect not 
unexpectedly becomes less and less pronounced with in- 
creasing perturbation, though the decay is always by or- 
ders of magnitude slower than for the case of a GOE 
perturbation (not shown). It is further seen that the lin- 
ear response approximation works very well up to about 
half the Heisenberg time, but underestimates the decay 
more and more for increasing r values. 

We thus can conclude that the fidelity freeze is not 
an artefact of the linear response approximation but is 
also present in the exact calculation. Due to the perfect 
and well established correspondence between random 
matrices and chaotic quantum systems ^|, [3 
this result provides an important new mechanism of 
preserving quantum stability. 
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APPENDIX A: CALCULATION OF 

Str [X{X + l)Fi] 

In terms of Pauli matrices X may be expressed as 



where 



X 



zl 



X 




where 



It follows 



X 



X + z 
y + z 



t1 + VGz 



(Al) 



(A2) 



Str[X(X + l)Pi] =4z(z-l) 



+ (l-2z)tr X(Pi),,. 



tr 



X2(P0,., (A3) 



where it was used that StrPi = StrP = 4, and where 
(Pi)^ y denotes the upper left submatrix of Pi. As was 
already mentioned, matrices Ui and U2 entering the cal- 
culation of Pi and P2 (see Eq. (|39|l ') are parameterized 
as 



Op 
1 



(A5) 



and Oi and O2 are 2x2 orthogonal matrices (see VWZ, 
Eq. (1.13)). The matrices Vp may be parametrized as 
(see VWZ, Eq. (K.26)) 



{Vpf' = 1 ± + \^'-^-'^Y'p ± Y^^-'^Y'p + 

(A6) 

where matrices Yi and Y2 are given by 



where 



Cp 



(A7) 



(A8) 



(see VWZ, Eqs. (K.23 + 25)). Note the convention 
(a*)* = —a for antisymmetric variables. In VWZ, Eq. 
(1.13) the sequence of the matrices on the right hand 
side of Eq. ljA4|l is reversed. Both parameterizations are 
equivalent and can be transformed into each other by a 
straightforward transformation of the ap, (3p variables. 

We are now going to calculate Pi = UiPUi^ — 
ViOiFOi^V^^ . To simplify notations, we shall omit 
the lower index '1' in the following. The calculation 
for P2 proceeds in the very same way. Since P = 
diag(l, 1,-1,-1) commutes with O, we are left with 



Pi = vpy-^ = V 



1 
1 



(A9) 



For the further calculation it is suitable to introduce 
the quantities 

A = aa* +P(3* , B = aa* - 13(3* , 

C = al3* +I3a* , D = i{a(3* - Pa*) . (AlO) 

A, B, C obey the relations 

A^^2a, B^ = -2a, C^ = -2a, = -2a , 

(All) 

where a — aa*(3P*, and 

AB AC = AD = BC = BD = CD = . (A12) 
It follows 

C^C - -Al - Ba, - Ca, , CC^ = Al (A13) 
As a direct consequence we have 



VpOp, (p=l,2) 



(A4) 



Y' 



yd 



-etc 
-CCt 

Al + Ba^ + Cax 
-Al 

-AY. 



(A14) 
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It follows from Eq. (|A6|I 
1 3 



V 



±1 



where 



w 



1 



8 

W ^LO* 
±U! W 



A\Y'± 1 



A 



A 
2" 



B 



C 



Y 



(A15) 



1 



,t - 



. A 3 



Pi 



{l-Aa + 2A)l + 2{Ba^ + Ca^) 2{l ~ A)C 



2(1 -A)C 



(-l-4a + 2A)l 



(A16) 



Inserting the results into Eq. ljA9p we have 



(A17) 



A corresponding expression is obtained for P2. In the 
average over the antisymmetric variables only the a terms 
survive, i. e. (Pi) = (P2) oc 1, as was stated above. 

Inserting finally the upper left corner element of Pi 
into Eq. HA3p . we end up with Eq. (|^ . 



APPENDIX B: CALCULATION OF 

Str[(l + 2X)Ki]V8 



It is suitable to write 
istr[(l + 2X)Ki]^ 



4str[(l + 2X),Ki]' 
Id 



1 



Str[(l + 2X)2Ki2] 



lstr[X,Ki]2 

+ istr[X(X + l)] 



(Bl) 



where Ki^ = = 1 was used. The second term on 
the right hand side is easily evaluated: 

^Stv[X[X + l)] = i[a;(x + l) + y(y+l) + 2z(l-z)] 



u(u + 1) + + z(l - z) 
-r^ + (2m + 1)t + . 



(B2) 



For the first term on the right hand side we need an 
expression for Ki. Using Ki = iJiKC/f ^ (see Eq. ^) 
and U = VO (see Eq. ^AA\ ) we may write 



(B3) 



since K (see equation ()20|l) commutes with O. Using Eq. 
HA15|I . we obtain 



Ki 



k 

K k 



(B4) 



where 

k = 



'WayW + uj'azOJ , 



U!(7yW — W(TzU! , 



(B5) 



Since Ki 
1 



1, we have 



^Str[X,Ki]^ = ^ Str(XKi)^-StrX2 



Str Xk 



StiX' 



, (B6) 



where in the second step expression (|A1|) for X was used. 
Using Eqs. (|A16|) and HB5|) we have 



(B7) 



fc= -{1 + A - D)(7y 



where C^^TzC = Da-y was used. Now the calculation of 
the terms entering the right hand side of Eq. I|B6() is 
straightforward: 



Istr(xfc)^ 



= (1 + 2A-2I?)(t 



2 v^) 



1 



StrX^ = 



(B8) 



Collecting the results of this subsection, we have 
^Str[(l + 2X)Ki]' = 

o 

t{2u + 1-t) -v^ + 2{A- D){t^ -v^) . (B9) 
whence follows Eq. 14711 . 
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FIG. 1: Ensemble average of the fidelity amplitude /e(r) with 
Ho taken from the GOE and a purely imaginary antisymmet- 
ric perturbation (solid line, calculated from Eq. H5()^ ') for 
different perturbation strengths e. For comparison the result 
from the linear response approximation (dashed line) , and for 
a GOE perturbation (dashed-dotted line) are shown as well. 
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FIG. 2: Ensemble average of the fidelity amplitude /c(r) for 
an imaginary antisymmetric perturbation as a function of e 
for three fixed values of t = 0.5, 1.0, 1.5 (from top to bot- 
tom, solid lines). Again the results from the linear response 
approximation are shown for comparison (dashed lines). 



